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ABSTRACT 


This  study  establishes  the  existence  of  damping  structu- 
ral properties  in  the  San  Andreas  Fault  zone  through  an  input 
model  of  surface  loading.   Any  changes  in  the  density,  mass 
or  bulk  elastic  properties  of  a  fault  zone  will  not  only  pro- 
duce a  change  in  the  gravity  field  of  the  earth,  but  also  a 
change  in  the  earth's  response  to  the  tidal  field  generated 
by  the  moon.   Knowledge  of  these  changes  is  important  to  en- 
gineering constructions,  pre-programmed  defense  systems,  etc. 
To  adequately  measure  these  responses,  sites  on  opposite 
sides  of  the  San  Andreas  Fault  zone  were  selected  and  tidal 
gravity  was  continuously  measured  for  finite  periods  of  time. 
The  Pacific  Ocean  and  its  very  low  frequency  lunar  portion  of 
the  ocean  load  spectrum  was  used  as  a  driving  function  to 
compare  observed  data  with  theoretically  generated  data  for 
the  same  sites.   A  least  squares  criteria  was  used  to  fit  the 
observed  and  theoretical  data  to  the  nine  most  important 
spectral  lines  and  to  obtain  the  phase  and  amplitude  for  each 
line  at  both  sites.   Utilizing  the  most  stable  lines  to  com- 
pare values,  a  significant  drop  (about  10$)  in  both  phase 
and  amplitude  occurred  between  the  two  sites.   The  anomalous 
decrease  in  both  areas  of  analysis  discloses  a  structural 
damping  in  the  San  Andreas  Fault  zone. 
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I.   INTRODUCTION 

A.   STATEMENT  OF  THE  PROBLEM 

The  combined  gravitational  fields  of  the  moon  and  the 
sun  acting  on  the  solid  earth  are  sufficient  to  displace  its 
surface  on  the  order  of  20  cm.  semi-daily  [l].   Unfortunate- 
ly, this  direct  forcing  function  is  so  broad  in  its  areal 
extent  (i.e.,  several  thousand  kilometers)  that  it  cannot 
reasonably  be  used  to  study  responses  of  small  crustal 
structures  (such  as  fault  zones)  [2]  which  are  on  the  order 
of  a  few  kilometers  in  width  and  depth.   However,  the  response 
of  the  ocean  to  the  lunar-solar  gravitational  field  (on  the 
order  of  meters)  constitutes  a  sufficiently  energetic  forc- 
ing function  such  that  under  certain  conditions  dynamic 
properties  of  crustal  structures  can  be  studied. 

The  conditions  under  which  the  ocean  can  effectively  be 
used  as  a  forcing  function  to  study  the  properties  of  a 
shallow  structure  are  defined  as: 

1.  The  distance  between  the  edge  of  the  ocean  load  and 
the  fault  structure  should  not  exceed  the  smallest  dimension 
of  the  fault  zone  by  one  order  of  magnitude  [3]; 

2.  Noise  response  in  the  earth  tide  records  obtained 
during  this  experiment  must  be  kept  to  a  minimum.   Minimiz- 
ing noise  response  is  achieved  by  using  parts  of  the  ocean 
tidal  spectrum  which  are  both  stationary  and  relatively  un- 
affected by  meteorological  influences. 


8 


The  two  lines  chosen  for  amplitude  and  phase  analysis  in 
this  study  are  the  principal  lunar  diurnal  0-,  and  the  princi- 
pal lunar  semi-diurnal  Mp  lines  which  are  the  most  energetic 
in  the  theoretical  tidal  spectrum  and  are  the  most  removed 
from  possible  contamination  due  to  meteorologically  induced 
noises. 

Earthquake  prediction  research  explores  the  causes  of 
earthquakes  from  two  points  of  view:   (l)  failure  brought 
about  by  an  external  increase  of  shear  stress  which  exceeds 
the  shear  modulus  of  the  fault;  (2)  failure  brought  about 
when  the  internal  properties  of  the  fault  change  (coefficient 
of  shear)  until  they  are  reduced  below  that  of  the  external 
shear  being  applied. 

In  this  study  responses  to  tidal  gravity  were  measured 
at  two  points  on  either  side  of  the  San  Andreas  Fault.   Tidal 
gravity  is  a  vertical  measure  of  the  acceleration  due  to 
gravity.   The  value  of  gravity  is  on  the  order  of  10  cm/sec/- 
sec  or  10  gal.   Resolution  on  the  order  of  10   gal  is 
achieved  with  the  instruments  used  (described  below)  giving 

Q 

results  accurate  to  one  part  in  a  10  .   The  distance  between 
the  two  stations  (about  5  kilometers)  is  small  compared  to 
the  distance  to  the  edge  of  the  ocean  load  (Pacific  Ocean, 
about  30  km).   The  amplitude  and  phase  of  the  two  measured 
lunar  tidal  constituents  (0-,  and  JVL)  were  compared  with  the 
amplitude  and  phase  of  the  same  two  lines  derived  from  an 
analysis  of  the  theoretical  tidal  gravity  calculated  at  each 


of  the  two  sites  for  identical  periods  of  time.   The  differ- 
ence in  the  observed  phases  between  the  sites  was  compared 
with  the  difference  in  theoretical  phases  between  the  sites 
and  a  10$  change  in  the  means  was  observed. 

This  result  suggests  that  the  principal  source  for  the 
phase  difference  between  the  observed  and  theoretical  values 
at  both  sites  is  due  to  the  response  of  these  sites  to  ocean 
tidal  loading  at  the  same  frequency.   However,  a  10$  devia- 
tion of  the  difference  between  these  two  sites  from  the  mean 
is  significantly  large  and  is  not  due  to  a  simple  shift  of  a 
few  kilometers  further  away  from  the  coastline.   Instead,  the 
anomalously  large  difference  in  the  comparison  is  explained 
by  the  fact  that  the  stations  are  not  just  separated  by  a 
few  kilometers  but  that  they  lie  on  opposite  sides  of  a  weak- 
ness in  the  earth's  crust.   The  absolute  difference  between 
the  most  inland  station  and  the  other  is  a  lag  of  about 
three  degrees  at  the  MU  frequency.   This  is  interpreted  pri- 
marily as  a  viscoelastic  response  of  the  fault  zone  to  ocean 
tidal  loading.   The  inland  edge  of  the  San  Andreas  Fault 
reaches  it  peak  amplitude  of  ocean  tidal  response  approxi- 
mately 1/10  hour  after  the  peak  response  of  the  ocean  side 
of  the  fault. 

While  this  study  did  not  attempt  to  directly  calculate 
any  geological  strength  parameters,  it  did  suggest  the 
existence  of  a  crustal  weakness  with  viscoelastic  properties 
that  strongly  dampen  the  M^  line  of  the  ocean  tidal  loading. 
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B.   BACKGROUND 

There  exists  in  any  multi-body  system  a  measurable  gravi- 
tational potential  due  to  the  attacting  fields  of  the  many 
masses  (or  bodies)  present.   Since  all  bodies  have  motion 
relative  to  all  others  in  a  system,  the  potential  fields 
fluctuate  with  time.   On  the  earth,  only  the  sun  and  the  moon 
contain  sufficiently  large  quantities  of  mass  at  a  close 
enough  distance  to  significantly  affect  tidal  variations  or 
time  fluctuations  in  the  potential  field. 

The  effects  from  the  sun  and  the  moon  generate  earth 
tides  in  a  similar  manner  to  their  production  of  ocean  tides. 
Unfortunately,  the  study  of  earth  tides  is  complicated  by 
the  presence  of  oceans  which  also  cause  changes  in  earth 
tides.   Although  quite  small,  ocean  loading  effects  on  the 
earth  tides  have  been  detected  in  measurements  as  far  inland 
as  mid-continent  in  the  United  States  [4]«   The  mass  of  the 
oceans  complicates  measurements  of  the  earth  tides  in  the 
same  ways  as  external  masses  (the  sun  and  moon)  plus   (1) 
the  loading  effect  due  to  water  (additional  mass)  and   (2) 
deformation  of  the  earth  due  to  the  redistribution  of  ocean 
waters  as 'the  water  more  freely  fluctuates. 

The  combination  of  effects  of  the  sun,  moon  and  oceans 
over  time  results  in  a  continuous  change  in  the  earth's 
gravity.   The  gravitational  force  at  any  given  point  in  time 
and  space  can  be  broken  into  two  components,  one  normal  to 
the  plane  tangent  to  the  earth's  surface  and  one  in  the 
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plane  tangent  to  the  surface.   This  study  will  concentrate 
on  the  former  component  called  tidal  gravity. 

Past  geological  and  seismological  studies  have  examined 
the  bulk  elastic  properties  of  the  whole  earth  and  the  dis- 
tribution of  these  properties  in  the  core  and  the  mantle. 
Few  studies  have  narrowed  their  objective  to  examining  the 
bulk  elastic  properties  of  only  the  crust  of  the  earth  and 
in  a  relatively  small  area.   Seismic  studies  show  that  the 
San  Andreas  Fault  zone  has  brittle  properties  to  a  depth  of 
only  ten  kilometers,  which  is  neither  the  bottom  of  the 
fault  nor  the  bottom  of  the  continental  plate  [5].   It  is 
approximately  one-third  of  the  depth  of  the  plate.   Authors 
of  these  studies  quite  often  generalize  their  results  to 
large  area  averages  and  tend  to  overshadow  the  relatively 
small  area  of  a  shallow  fault  zone  in  the  final  analysis  [2], 
The  seismic  studies  and  their  related  high  frequencies  (from 
earthquakes)  also  do  not  allow  determination  of  anelastic 
responses,  but  very  low  frequency  studies  such  as  those  re- 
sulting from  ocean  loading  should  disclose  any  anelastic 
properties  present. 

This  study  is  a  portion  of  a  larger  project,  (headed  by 
Dr.  M.  D.  Wood)  at  the  National  Center  for  Earthquake  Re- 
search (NCER),  to  establish  that  the  response  of  the  earth 
to  tidal  loading  is  significantly  perturbed  by  the  San 
Andreas  Fault  and  also  to  determine  through  spatial  and 
temporal  means  the  bulk  properties  of  the  fault  zone. 
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There  are  two  possible  sources  of  error  in  analyzing 
these  data  that  should  be  noted.   First,  due  to  many  con- 
straints, the  three  experiments  in  this  study  were  performed 
sequentially  instead  of  simultaneously.   Second,  the  San 
Francisco  Bay,  due  to  its  size  and  its  proximity  to  the  sites, 
will  affect  the  values  measured  in  the  experiment.   The  ef- 
fects of  these  errors  are  assumed  to  be  relatively  small, 
but  their  exact  effects  are  undetermined.  . 
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II.   DATA  ACQUISITION  AND  PROCESSING 

A.   LOCATION  SELECTION 

The  locations  selected  for  the  two  main  sites  for  this 
analysis  were  obtained  from  a  threefold  criterion.   First, 
it  was  desired  to  ensure  that  the  two  sites  were  on  opposite 
sides  of  the  San  Andreas  Fault  zone  and  nearly  equidistant 
from  the  fault  region.   This  was  necessary  in  order  to  ob- 
tain spatial  analysis  of  the  tidal  wave  as  it  passed  from  one 
site  to  the  other  by  traversing  through  the  fault  region. 
Second,  the  sites  were  desired  to  be  on  a  line  nearly  per- 
pendicular to  the  coastline.   Third  was  the  availability  of 
suitable  sites.   It  was  imperative  that  the  physical  location 
of  each  instrument  be  based  on  a  very  sturdy  foundation, 
ideally  on  deep  lying  rock,  to  minimize  noise.   To  further 
eliminate  noises  (in  the  form  of  other  energies  such  as  wind, 
temperature  fluctuation,  surface  vibrations,  etc.)  it  was 
desirable  to  place  each  recording  instrument  below  the  sur- 
face of  the  earth.   Time  and  financial  restrictions  forced  a 
compromise  of  this  desire  to  a  level,  surface  (or  very 
shallow)  base  that  was  protected  from  as  many  of  the  other 
noise-producing  sources  as  possible.   In  addition  to  these 
restrictions,  the  security  of  the  expensive  recording  instru- 
ments became  a  major  concern.   The  instruments  could  not  be 
left  in  the  open  where  vandals  or  sightseers  might  tamper 
with  them.   It  was  also  desired  to  have  daily  time  marks 
entered  on  each  record,  generating  the  requirement  for 
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someone  to  be  present  at  the  site  each  day  to  perform  this 
task.   Appropriate  sites  with  the  necessary  attendance  re- 
quirements were  located  [Figure  1].   One  recording  instru- 
ment was  in  a  basement  vault  on  a  base  on  concrete  poured  on 
solid  rock  (Block  experiment);  another  was  in  a  separate 
garage  on  a  concrete  base  poured  on  rock  (Rose  experiment). 
Both  properties  were  private  and  actual  instrument  locations 
were  secure.   The  third  site  was  located  in  the  NCER  Build- 
ing in  Menlo  Park  on  a  concrete  floor  (NCER  experiment). 

B.  EARTH    TIDE  RECORDINGS 

All  the  earth  tide  recordings  were  taken  from  LaCoste- 
Romberg  Earth  Tide  Gravity  Meter  type  instruments  [l]  on 
paper  strip  charts  [Figure  2],   These  were  recorded  from 
NCER  owned  instruments  run  by  120v  AC  with  I2v  automobile 
batteries  as  a  back-up. 

To  allow  future  adjustments  for  inaccuracies  in  chart 
speed,  stretching  or  shrinking,  etc.,  daily  time  marks  were 
mechanically  entered  on  each  chart  and  labeled.   Chart  speed 
for  each  experiment  was  two  inches  per  hour  with  a  displace- 
ment of  one  inch  equal  to  30  /u,  gal.   Each  experiment  was 
attempted  to  run  at  least  thirty  days  to  give  a  large  enough 
sample  for  adequate  analysis  [6]. 

C.  DATA  PROCESSING 

The  strip  charts  were  digitized  onto  magnetic  tape  by 
use  of  a  Calma  302  digitizer.   The  cross-hair  cursor  moved 
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Figure  1. I  Test  Site  Locations 
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vertically  or  horizontally  as  desired  across  the  strip 
chart  taking  discrete  samples  every  .01  inches  and  recording 
the  readings.   A  special  indication  was  noted  at  each  daily 
time  mark  entered  on  the  chart  by  monitoring  personnel,  and 
the  time  and  date  were  recorded  by  the  person  performing  the 
digitization.   The  digital  tape  data  was  later  corrected 
for  chart  stretching,  shrinkage,  etc.,  by  proportionally 
distributing  the  actual  recorded  values  between  time  marks 
to  the  correct  number  of  .01  inch  increments  (.01  inch 
corresponds  to  IS  seconds  of  time)  required  between  the  daily 
time  marks.   This  localized  physical  errors  between  time 
marks  (usually  every  24  hours)  and  minimized  accumulation 
of  error. 

The  digitized  data  consisted  of  readings  based  on  a 
continuous  time  record.   These  data  were  then  processed  with 
a  computer  program  (designed  by  M.  D.  Wood  and  P.  A.  Marshall 
[7]  from  the  National  Center  for  Earthquake  Research  (NCER) 
and  modified  by  M.  D.  Wood  and  R.  L.  Keller  for  use  in  this 
study)  which  decimated  the  sample  interval  from  .01  inches 
to  hourly  values  taken  on  the  hour. 

Sampling  theory  requires  a  sample  interval  at  least  twice 
as  often  as  the  highest  expected  frequency.   Implicit  in  this 
is  the  assumption  that  the  recorded  data  be  of  a  reasonable 
duration  of  time.   A  reasonable  time  duration  is  propor- 
tional to  the  reciprocal  of  the  smallest  difference  between 
any  two  frequencies  expected  in  the  data  [6].   Since  only 
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nine  lines  of  the  SO-plus  theoretical  lines  are  used,  it  can 
be  calculated  that  at  least  one  month  of  data  is  necessary 
for  a  Vfo   resolution  of  tidal  data  [6],   The  smallest  fre- 
quency line  used  in  this  study  is  the  M,  with  a  frequency 
of  less  than  4  cycles/day,  requiring  at  least  8   samples  per 
day.   By  oversampling,  the  resolution  remains  well  within  Vfo. 

III.   LEAST  SQUARES  ANALYSIS 

A.   AMPLITUDE  ANALYSIS 

The  observed  tidal  gravity  records  were  fit  with  a  least 
squares  computer  program  to  the  nine  theoretical  lines  of  the 
earth  tide  spectrum  (explained  in  Appendix  A).  The  computer 
program  (designed  by  M.  D.  Wood  [6]  and  modified  by  M.  D. 

Wood  and  N.  E.  King  of  NCER  and  R.  L.  Keller)  performs  a  re- 

th 
gression  of  the  data  to  fit  both  an  n —  polynomial  and  a 

sine  wave  function  simultaneously.   The  results  give  the 
amplitudes  and  phases  of  the  nine  lines  mentioned  above. 

It  is  a  well  established  fact  that  the  principal  lunar 
semi-diurnal  IVL  and  the  principal  lunar  diurnal  0-.  lines  of 
the  earth  tide  frequency  spectrum  are  the  most  invariant  of 
all  the  tidal  lines  [l],  [6],  [£].   Of  these  two,  the  Mg 
line  has  proven  to  be  less  variant  from  its  theoretical  value 
than  the  other  lines  in  previous  studies  conducted  in  this 
region  [6],  [7]»  [&]  and  was  thus  used  as  a  normalizing  fac- 
tor for  all  of  the  spectral  lines.   The  normalized  amplitudes 
were  computed  using  the  following  formula: 
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Ai  =  Ai/AM2 

where  A.  is  the  (dimensionless )  relative  amplitude,  A.  is 
the  amplitude  of  the  constituent  line  being  normalized,  and 
AM2  is  the  amplitude  of  the  M^  line.   This  normalizing  pro- 
cess results  in  dimensionless  values,  facilitating  compari- 
son between  observed  and  theoretical  results. 

The  normalizing  of  the  amplitudes  of  the  selected 
lines  in  the  earth  tide  spectrum  was  conducted  for  both  the 
theoretical  data  [Table  III]  and  for  the  observed  data  [Table 
IV].   Visual  displays  of  these  data  are  given  in  Figure  7 
and  Figure  8. 

Taking  the  expected  amplitude  (theoretical)  of  each 
constituent  line  for  each  of  the  sites  and  subtracting  the 
corresponding  observed  amplitude,  the  amplitude  differences 
are  obtained  [Table  I].   Differences  can  be  expected  due  to 
forcing  nine  lines  of  a  frequency  spectrum  to  absorb  what  is 
acknowledged  to  consist  of  at  least  $9  lines  and  also  due  to 
other  non-tidal  factors  such  as  seasonal  noise,  extreme  wea- 
ther conditions,  etc.   These  differences  should  remain 
about  constant  between  sites  (resulting  in  low  values  in  the 
last  column)  and  would  indicate  a  high  stability  of  a  con- 
stituent line.   The  last  column  of  Table  I  shows  a  standard 
deviation  of  about  5^  for  the  0-,  line  from  its  average  theo- 
retical mean.   This  indicates  that  the  0-,  line  could  be  pre- 
dicted to  have  a  measured  (observed)  value  with  a  standard 
deviation  equal  to  about  5$  of  the  theoretical  value  of  the 

amplitude. 
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The  normalizing  of  the  amplitudes  of  the  selected 
lines  in  the  earth  tide  spectrum  was  conducted  for  both  the 
theoretical  data  [Table  II]  and  for  the  observed  data  [Table 
III].   Visual  displays  of  these  data  are  given  in  Figure  3 
and  Figure  4.   These  tables  and  figures  are  presented  here 
to  facilitate  comparison  of  the  results  obtained  in  this 
study  to  those  of  other  completed  studies  [6],  [7],  [&] 
making  use  of  similar  types  of  displays. 

A  relative  gravimetric  factor  G.  can  be  computed  by 
the  following  formula: 

G.   =  A.  ( experimental) /A .  (theoretical) 

where  the  A.  are  described  above.   Using  the  two  most  stable 

l 

lines  to  establish  a  comparative  base,  we  compute: 
Gr(01/M2)  =  Ar(01/M2)(  experimental )/Ar( 0-,/^) (theoretical) 

Theoretically,  this  ratio  should  equal  to  1.00,  but  with 

other  dynamic  considerations  a  Vfo   variation  [9]  is  accepted. 

Ten  world-wide  stations  studied  by  Harrison  and  others  [10] 
resulted  in  an  average  value  of  : 

Gr(01/M2)  =  0.9^9+0.035. 

Four  stations  in  the  San  Francisco  Bay  area  analyzed  by  Wood 
[6]  gave  an  average  value  of: 

Gr(01A2)  =  0.933+0.026. 
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The  average  value  obtained  by  the  three  stations   in  this 
study  gave: 


G x(01/k2)    =  1.051+0.061, 


These  values  should  not  be  different,  regardless  of  loca- 
tion, yet  this  study  and  that  by  Wood  indicate  values  about 
ki°   to  5i°   outside  of  the  accepted  range.   Both  of  these 
studies  were  conducted  on  the  West  Coast  near  the  San 
Francisco  Bay  Region  and  could  indicate  a  local  disparity 
from  areas  with  thicker  crustal  depths  (The  California  coast- 
al regions  are  known  to  have  a  much  thinner  crustal  depth 
than  most  other  continental  areas).   However,  no  conclusion 
is  drawn  from  this  study  other  than  to  mention  that  this  dis- 
parity should  not  exist. 

B.   PHASE  ANALYSIS 

The  observed  tidal  gravity  records  were  simultaneously 
fit  in  the  least  squares  program  (described  in  the  previous 
section)  to  the  nine  theoretical  waves  of  the  earth  tide 
spectrum  to  estimate  the  phase  of  each  constituent  line. 
The  values  for  both  theoretical  and  observed  phases  for  the 
IVU  lines  are  given  in  Table  IV.   All  charts  were  maintained 
in  reference  to  Greenwich  Time  (GMT)  except  for  the  Rose 
site  which  was  marked  in  local  daylight  time.   A  seven  hour 
difference  existed  between  these  base  reference  times  and 
the  phase  was  appropriately  advanced  for  this  discrepancy. 
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The  result:,  show  a  phase  difference  between  the  theoreti- 
cal and  observed  M«  line  at  both  the  Rose  r.ite  and  the  Block 
site. 
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Figure   3.      Composite  Theoretical  Amplitude  Histogram 
Normalized  by  NL   Component 
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Figure  4.   Composite  Observed  Amplitude  Histogram, 
Normalized  by  Mp  Component. 
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IV.   RESULTS 

The  average  value  of  G  (0-,/Mp)  was  much  higher  than  those 
produced  by  the  experiments  done  by  Harrison  and  others  [10] 
and  Wood  [7].   As  noted  above,  no  explanation  for  this  has 
to  date  been  found.   The  relative  gravimetric  factor  showed 
a  downward  change  between  the  two  sites.   It  decreased  by 
more  than  0.1,  or  about  a  10$  drop  in  value. 

In  the  examination  of  the  two  most  stable  lines  in  the 
earth  tide  spectrum,  the  value  of  the  0-,/Mp  (the  normalized 
amplitude  of  the  0-,  constituent  line)  decreased  from  the  Rose 
to  the  Block  station  by  a  value  of  .04$$  which  is  about  a 
7M°   decrease  between  the  two  sites.   Theoretically,  the  nor- 
malized 0-,  line  amplitude  should  not  show  much  change  in 
value,  specifically  across  a  distance  of  5  km.   From  the 
values  in  Table  II,  it  can  be  shown  that  about  a  2fo   increase 
should  have  occurred. 

There  also  exists  a  phase  difference  between  the  theore- 
tical and  observed  data  at  each  site.   The  observed  and 
theoretical  lines  should  all  be  in  phase  with  each  other. 
Allowing  for  small  variations  due  to  noise  and  tectonic  in- 
terferences prior  to  the  arrival  of  the  M^  line  at  the  Rose 
site  (again  the  most  stable  line  is  selected  for  comparison), 
there  may  exist  a  phase  difference  between  theoretical  and 
observed  lines.   Between  the  Rose  and  Block  sites,  however, 
this  difference  should  remain  constant.   In  fact,  it  does 
not;  it  changes  by  more  than  2-$   degrees.   This  change  is  a 
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decrease  in  phase  difference,  or  a  lag,  of  about  9"2^  between 
the  two  sites.   Physically  interpreted,  this  is  a  time  delay 
of  about  1/10  hour  for  peak  amplitude  of  the  ]VU  wave  to 
travel  from  the  Rose  to  the  Block  site. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

With  the  evidence  advanced  in  this  study,  it  appears 
that  some  disturbing  factor  has  caused  a  large  percentage 
decrease  in  tidal  amplitude  and  phase  between  the  Rose  and 
the  Block  sites.   These  percentage  changes  are  much  too  great 
to  be  explained  by  the  5  km.  distance  between  the  sites.   The 
fact  that  the  San  Andreas  Fault  is  located  between  these  two 
sites  leads  one  to  suspect  that  drastic  changes  in  earth 
tides  between  the  two  sites  might  be  due  to  the  fault  region 
and  its  different  structural  properties  from  surrounding  area 
properties. 

Assuming  that  the  fault  region  is  responsible  for  these 
tidal  changes,  it  is  possible  to  calculate  the  approximate 
strength  parameters  of  the  geological  structure  of  the  fault 
region  itself  through  a  more  detailed  geological  study  of 
sites  across  the  fault. 

Since  this  experiment  was  limited  in  time,  money,  equip- 
ment, available  sites,  etc.,  it  cannot  provide  absolute 
proof  that  the  fault  region  is  responsible  for  the  observed 
changes.   The  study  provides  evidence  that  the  fault  might  be 
reason  for  the  observed  changes.   It  is  recommended  that 
further  investigations  of  a  similar  nature  be  conducted  at 
sites  void  of  an  inland  bay  to  verify  the  results  of  this 
study.   These  future  studies  could  also  be  set  up  to  simul- 
taneously examine  the  specific  geological  properties  exhibited 
by  the  fault  region  in  each  area  of  study. 
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APPENDIX  A 
THEORY  OF  EARTH  TIDES 

A  study  of  the  earth  tide  curve  shown  in  Figure  5  reveals 
the  presence  of  at  least  two  frequency  components  of  the 
earth  tide  spectrum.   Identified  in  this  figure  are  the  diur- 
nal frequency  with  a  period  of  approximately  2$  days  and  the 
semi-diurnal  frequency  with  an  approximate  14  day  period.   A 
Fourier  analysis  of  a  record  such  as  this  reveals  that  there 
exist  many  constituent  waves  that  combine  to  give  us  the  wave 
form.   Analysis  of  longer  records  would  reveal  semi-annual, 
annual,  and  even  longer  period  waves  which  correspond  to  mo- 
tions of  the  sun  and  moon  in  relation  to  the  earth. 

When  the  complex  relationships  between  the  moon  (and  sun) 
of  a  rigid  earth  are  considered,  the  constituent  waves  (lines 
in  the  spectrum)  derived  via  the  Fourier  analysis  become 
meaningful.   Let  us  examine  a  fictitious  model  [6],  [$]  in 
which  the  moon  is  orbiting  about  the  earth  in  the  plane  of 
the  earth's  equator  at  a  constant  distance  with  a  constant 
velocity  equal  to  the  mean  distance  and  velocity  of  the  actu- 
al moon.   The  earth  is  also  revolving  on  its  axis  once  every 
24  hours  and  this  motion  plus  the  moon's  gravitational  at- 
traction on  the  earth's  surface  produces  a  tide  with  a  per- 
iod equal  to  a  lunar  half-day  and  is  described  by  the  IVL 
principal  lunar  semi-diurnal  line,  the  most  energetic  of  the 
semi-diurnal  lines.   In  similar  fashion,  the  S2  principal 
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solar  semi-diurnal  line  is  produced  by  a  fictitious  sun  in  a 
circular  orbit  at  a  constant  velocity  and  distance. 

If  the  moon's  declination  in  its  orbit  (above  or  below 
the  plane  of  the  earth's  equator)  is  now  added  to  our  simple 
model,  the  distance  between  the  moon  and  an  observation  point 
on  the  earth  varies  with  time,  but  in  a  different  manner  than 
the  simpler  model.   In  the  simple  model,  the  distance  varia- 
tion with  time  was  due  only  to  the  observer's  position  and 
the  lunar  position  in  the  equatorial  plane.   Now  a  declina- 
tion distance  away  from  the  equatorial  plane  also  varies  with 
time  as  well  as  the  consideration  for  the  earlier  model,  add- 
ing another  variable  to  be  considered  in  the  distance  varia- 
tion between  the  moon  and  an  observation  point  on  earth. 
For  observation  points  on  the  equator,  the  effects  of  the 
lunar  orbit  in  the  northern  and  southern  hemisphere  are  equal 
in  magnitude.   For  non-equatorial  observers,  the  amplitude 
of  the  two  diurnal  maxima  would  vary  in  time  with  the  chang- 
ing distance  from  the  observation  point  to  the  moon.   This 
is  referred  to  as  the  diurnal  inequality.   The  cyclical  mo- 
tion of  the  moon's  declination  causes  a  fortnightly  modula- 
tion of  the  earth  tides. 

Adding  the  ellipticity  of  the  moon's  orbit  to  our  model 
now  leads  to  the  cyclical  changes  in  the  distance  from  the 
moon  to  the  earth  and  causes  modulation  of  both  the  amplitudes 
and  phases  of  all  the  above  effects. 
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Associating  a  theoretical  wave  with  each  of  the  motions 
of  the  moon  does  not  complete  an  explanation  of  the  earth 
tide  spectrum.   Secondary  waves  are  apparent  due  to  the  beat 
frequencies  of  the  principal  tide  waves.   These  beat  fre- 
quencies are  obtained  from  the  sum  and  difference  frequencies 
of  the  two  principal  waves.   When  the  difference  in  the  lunar 
and  solar  principal  semi-diurnal  waves  (M^  and  Sp)  is  taken, 
a  fortnightly  modulation  occurs  at  the  beat  frequency.   The 
lunar  elliptic  semi-diurnal  wave,   No  >  and  the  JVU  wave  com- 
bine for  a  beat  frequency  to  form  a  modulation  with  a  monthly 
period  that  is  associated  with  the  eccentricity  of  the  lunar 
orbit.   The  K2  luni-solar  semi-diurnal  declinational  wave  is 
composed  of  lunar  and  solar  components  which  are  analytically 
inseparable,  as  is  the  K-,  luni-solar  diurnal  declination  wave. 

Further  beat  frequencies  exist  which  yield  waves  derived 
in  a  similar  fashion.   Table  I  lists  the  principal  theoreti- 
cal- tidal  waves  of  the  earth  tide  spectrum  which  are  of 
importance  to  this  study.   The  amplitude  of  each  of  these 
waves  is  dependent  on  latitude  and  varies  with  time. 
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APPENDIX  B 
THEORETICAL  TIDE  GENERATION 

The  gravitational  potential  at  some  point  P  on  the  sur- 
face of  the  earth  [Figure  6]  due  to  another  mass  (moon  or 

sun)  can  be  expressed  as  [6],  [11],  [12]: 

2  f^.  /-\-  oo 


w  =  m^y  y&J     •   Pn(z)  =    Y-     w 

n=o  n=o 


n 


where    G    is  a  gravitational  constant 
M    is  the  external  mass 
r    is  the  radius  of  the  earth 
R    is  the  absolute  distance  from  the  earth 

to  the  external  mass 
P  (z)  is  the  Legendre  polynomial 

z    is  the  angular  position  of  the  external 

mass  relative  to  the  zenith  of  point  P 
P    is  a  point  on  the  earth's  surface. 

Expanding  the  terms  in  the  summation  above  necessitates  the 

calculation  of  the  Legendre  polynomials.   The  first  few  terms 

of  the  Legendre  polynomials  are  [12]: 

P0(z)=l 

Px(z)  =  z 

P2(z)  =  i(3z2-l) 

P3(z)  =  i(5z3-3  ) 

P4(z)  =  l/8(35z4-30z2+3) 

Examining  these  results  it  is  seen  that  the  first  term, 
W  =GM(l/R),  is  constant  and  has  no  effect  on  changing  the  gravi- 
tational  potential.   The  second  term,  W-,  =GM(  r/R  )cos(z),  uni- 
formly affects  all  points  on  the  surface  of  the  earth  and  has 
no  effect  on  changing  the  potential.   The  third  term  (and  all 
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remaining  terms)  affect  the  tidal  potential  and  are  all 

candidates  for  significant  effects  in  the  summation.   The 

r/R  term  in  the  summation  is  approximately  1/60  for  the  moon 

and  1/2400  for  the  sun  and  after  n=3  accounts  for  less  than 

2$  of  the  previous  term.   Thus,  the  gravitational  potential 

can  be  reduced  to  the  simple  expression: 
2  2 

W2  =  GM^I  P2(z)  =  GM^J  !-(3cos2(z)-l). 

To  express  z  in  terms  of  fundamental  angular  parameters 
of  position,  spherical  trigonometry  is  applied  where: 

6   =  declination  of  external  mass  M 

t  =  hour  angle  of  point  P 

<*  =  zenith  angle  of  external  mass  at  point  P 
f  =  latitude  of  P 
X     =   longitude  of  P 

z  =  polar  angle  of  external  mass  M 
Using  the  trigonometric  relation: 

cos(cc)=   cos(e)cos(z)  +  sin(e)sin(  z)cos(t ) 
and      cos(&)=  sin(?  ),  sin(e)  =  cos  (?) 

cos(z  )=  sin(  6  ),  sin(z  )  =  cos  (S  ) 

the  gravitational  potential  becomes 

2  ? 

W2=GM£j  i  (3[sin(  f  )  sin(  6  )  +  cos(f  )cos( 6  )cos(t)]  -1) 
R 

and  simplifies  to 

W2  =  3K(sin2(f)-l/3)(sin2(6)-l/3) 
+  K(sin(2f)sin(2S)cos(t)) 

+  K(cos2(£)cos2(6)cos(2t)) 


where  K  =  3/4GMr 


This  equation  is  shown  in  Figure  7. 

Using  the  moon  as  an  example,  the  first  term  in  this  ex- 
pression contains  the  variable  6  ( f   is  constant  for  a  given 
point  P) ,  namely  sin  (6)  ,  which  is  an  even  function  varying 
as  the  declination  of  the  mass  M  traverses  from  its  maximum 
positive  value  to  its  maximum  negative  value  and  back  each 
month.   The  functional  value  stays  positive  (since  sin  (<5  ) 
is  an  even  function)  and  the  frequency  of  the  first  term  be- 
comes two  cycles  per  month  and  is  called  the  fortnightly  wave, 
The  second  term  varies  with  6  and  t,  but  8  completes  a  cycle 
in  one  month  while  t  has  a  period  of  about  one  day.   There- 
fore, this  term  has  a  frequency  of  one  cycle  per  lunar  day 
and  is  called  the  diurnal  wave.   Similarly,  the  third  term, 
with  a  frequency  of  two  cycles  per  lunar  day,  is  called  the 
semi-diurnal  wave. 

A  similar  situation  occurs  with  the  sun  as  the  external 
mass,  and  the  sum  of  the  results  from  the  moon  and  the  sun 
yields  the  terms  of  the  tidal  potential.   The  remainder  of 
the  tidal  potential  spectrum  consists  of  the  higher  ordered 
terms  and  combined  lunar-solar  terms  neglected  in  this  analy- 
sis.  For  the  purposes  of  this  study  the  nine  waves  described 
in  Table  V  are  sufficient  since  they  contain  more  than  9&/0 
of  the  effects  of  the  39  theoretical  lines. 

In  the  generation  of  the  theoretical  earth  tides  for 
this  study,  J.  C.  Harrison's  program  [13]  was  used.   In  the 
development  of  his  model  for  generating  earth  tides,  Harrison 
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uses  Wp  for  the  moon  (as  described  above)  and  Wp  for  the  sun 
and  also  includes  W^  for  the  moon  to  achieve  greater  preci- 


sion in  his  results. 
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Figure  6.   Geometric  Relation  between  Station  Coordinates  and 
Tide  Producing  Body  Angular  Coordinates.   Where  <*• 
is  the  zenith  angle  of  the  tide  producing  body  at 
the  station  location  P.   P  is  at  colatitude $  and 
longitude  X .   Z  is  the  polar  angle  and  f  the 
terrestrial  east  longitude  of  M.   (Adopted  from 
Munk  and  Cartwright) 
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